Disordered Carbon nanotube alloys in the Effect Medium Super Cell Approximation. 
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We investigate a disordered single-walled carbon nanotube (SWCNT) in an effective medium su- 
per cell approximation (EMSCA). First type of disorder that we consider is the presence of vacancies. 
Our results show that the vacancies induce some bound states on their neighbor host sites, leading 
to the creation of a band around the Fermi energy in the SWCNT average density of states. Second 
type of disorder considered is a substitutional B c bN cn Ci- c b-cn alloy due to it's applications in het- 
rojunctions. We found that for a fixed boron (nitrogen) concentration, by increasing the nitrogen 
(boron) concentration the averaged semiconducting gap, E g , decreases and at a critical concentra- 
tion it disappears. A consequence of our results for nano electronic devices is that by changing 
the boron(nitrogen) concentration, one can make a semiconductor SWCNT with a pre-determined 
energy gap. 



The role of disorder in a SWCNT is of importance from 
two prospectives; first in the growth process of a SWCNT 
due to the experimental environment some impurity 
atoms are inserted and vacancies are created 0, U 0- 
Second, we deliberately implant the impurity so as to 
construct a new nanotube alloys, such as BrhC\- c b, 
Ncnd-cn and B cb N cn Ci- cb - cn SWCNTsHSEl, with 
pre-determined physical properties. In the first case, the 
effect of a point-like defect was investigated by calcu- 
lation of electron reflection coefficient p], and also two 
substitute defects in an armchair SWCNT Q .It has been 
found that the symmetry of defects strongly affected 
the conductance and the local density of states. By 
different techniques, the boron nitride SWCNTs junc- 
tions 0, H(3 | , the spin polarization in a quasi one dimen- 
sional C/BN nanotube[ll| and also the current distribu- 
tion in boron and nitrogen doped SWCNTs were inves- 
tigated |l2^ . For finite impurity concentration, a system- 
atic field theory technique beyond single-site T-matrix 
approximation has not yet been applied to the disor- 
dered SWCNTs [13 - In this Paper, for the first time, 
by applying the EMSCA0 El] method to the disor- 
dered SWCNT, we will go beyond this approximation 
and consider the presence of finite impurities. We pro- 
vide a more realistic description of the effects of disorder, 
due to vacancies, on an armchair SWCNT's and a zigzag 
SWCNT's density of states (DOS). Also in this formal- 
ism, we address the question of how the doping of a zigzag 
SWCNT by boron (nitrogen), i.e. B cb N cn Ci- c b- C m con- 
trols the semiconducting gap, E g . 

Let us consider the Hamiltonian as a general random 
tight-binding model0] 



H 



where t"f a<7 are the hopping integrals between the tt or- 



£ 



t aP c c 



t 

■ c ■ 

1(7 j a 



(i) 



ijaa 



bitals of sites i and j with spin a. a and f3 refer to the 
A or B sites, [i is the chemical potential and ef is the 



random on-site energy where it takes with probability 
1 — c for host sites and S with probability c for impurity 
sites. For the B cb N cn Ci^ cb - cn SWCNT alloy, ef takes 5 
(— t) with probability cb for boron sites, -8 (= — t) with 
probability cn for nitrogen sites and with probability 
1 — cb — cn for the carbon sites Q, where t is hopping 
integral to the nearest neighbour. The matrix form of 
EqUJis, 

H = - £ '"-vv + E - ( 2 ) 
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where the two-component field operator, is given by 

= i ) ■ (3) 
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and £j is the random on-site energy matrix, 
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and iij is the hopping matrix defined by, 
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and I is a 2 x 2 unitary matrix. 

The equation of motion for electrons in such a lattice 

is, 

( (EI - §i + fii)5ii - hi ) G(Z, j; E) = IS l3 (6) 



where G(i,j;E) is the random Green function matrix 
defined by 

C(4 i-F\- ( GAA ^r,E) G AB (i,j;E) \ 
G{t,j,E)-\ G BA {iJ . E) G BB (iJ . E) ■ (7) 
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We considered the ii as a perturbation parameter, 
hence G(i,j;E) in Eq0 may be expanded in terms of 
the perfect Green function matrix G°(i,j;E) as, 

G(i,j; E) = G°(i, j; E) + Y G°(i, I; E)iiG(l,j; E) (8) 

i 

where G (i,j; E) is given by 

G°(i,j;E) = (EI- ek + I^)- 1 . (9) 

k 

with ek = ]7 Sij % eSk ' rM being the band structure for 
perfect system. In our calculations we assumed allowed 
hopping to the nearest neighbors and neglected the oth- 
ers. Hence 
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and the dispersion relation is 
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where 7(k) = X)i=i e jk Ti and are three vectors that 
connect an A(B) site to it's nearest neighbors B(A) sites. 

The Dyson equation for the averaged Green function, 
G(i,j;E), corresponding to Eq|S]is 

G(i,j;E) = G°(i,j; E) 

+ J2 G °(^ l ' E )*( l > l';E)G(l',j; E),(12) 
w 

where the self energy E(Z, Z ;E) is defined by 

(e l G(l,j;E))=Y,X(l,l';E)G(l\j;E). (13) 
i' 

The Fourier transform of G(i, j; E) in Eg 1121 is given by 
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S(k; J B) = -^e- lk ' r «S(i ! i;^), (15) 



is the self energy Fourier transform. 

We solve Eq|8] using the EMSCA method Hi EH for 
the case of four sites super cell, i.e. N c = 4. FigQ] shows 
a 2-dimensional graphene sheet. Each cell of the Bravias 
lattice includes two nonequivalent sites that are denoted 
by A and B. The primitive vectors of the Bravias lat- 
tice are a and b and the chiral vector is L. The heavy 
dashed line on the figure shows a four-sites super cell of 
the graphene Bravias lattice. 




FIG. 1: A two dimensional graphene sheet. The light dashed 
lines illustrate the Bravias lattice unit cells, a and b are the 
primitive vectors. Each cell includes two non-equivalent sites, 
which are denoted by A and B. L — n a a + nth is the chiral 
vector. For an armchair SWCNT n a — n and rib = 2n, while 
for a zigzag SWCNT n a = n and n& = 0. The heavy dashed 
line denotes a four-site super cell. 



In the EMSCA technique, the super cell random Green 
functions, G*™(i,j;E), are related to the cavity Green 
function Q(i,j;E) via 



(11) G£(I, J; E) = Q(I, J; E)+Y 0(1, L; E)e L G™(L, J; E), 



(16) 

where {7} refers to the sites inside the super cell. Also 
the Dyson's-like equation for the average super cell Green 
function, G SC (I, J; E), is given by 

G SC (I,J;E) = g(I,J;E) 

+ J20(I,L;E)V SC (L,L;E)G(L,J;E). 
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The Fourier transform of G sc (7, J; E) in Eall7lis 



G SC (K„; E) = g(K n ; E)+g(K n ; E)V sc (K n ; E)G(K n ; E), 

(18) 

where 

E(K n ; E) = -i- e K "' riJ S(7, J; E) (19) 
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and 

G(K n ;E) = ^ 

x Y, (lS-e k + I/x-S(K n ;£)) _1 . 

k£ nth patches 

(20) 

To calculate the G SC (I, J; E) and G™ P (I, J; E), Eqsd 
[30] should be solved self consistently. 

A SWCNT with vacancies is considered, the aver- 
aged density of states for different vacancy concentra- 
tions is calculated. We found that vacancies create some 
bound states around the Fermi level on their host neigh- 
bour sites, hence constructing a band in the averaged 
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FIG. 2: (a) Comparison of average density of states for 
a: (10,10) armchair SWCNT and (b): for a (10,0) zigzag 
SWCNT at S = 40t and different impurity concentration. At 
high random energies, S >> band width, we have band split- 
ting and the impurity band is located at higher energies, while 
remaining sites are the vacancies that induce a band around 
the Fermi energy. By increasing the vacancy concentration 
the height of the peak at the Fermi energy increases and the 
van How singularities are smeared out. 

density of states. Also the one-dimensional (ID) van 
How singularities in high vacancy concentrations disap- 
pear. FigJ2fa),(b) shows the comparison between the av- 
erage density of states for different vacancy concentra- 
tions in (10,10) and (10,0) SWCNTs respectively. The 
bound states due to vacancies around the Fermi energy is 
marked by an arrow. In short, our results show that va- 
cancies not only change the average density of states but 
also the number of electrons located on the host sites and 
also at high vacancy concentrations SWCNT's loses their 
ID characteristics and become similar to a 2D disordered 
geraphene sheet. 

We now investigate the effect of nitrogen and boron 
doping on a (10,0) zigzag SWCNT. Two cases are con- 
sidered, first fixed boron concentration at cb = 0.15, with 
variable nitrogen concentration. In this case, we found 
that the average semiconducting gap, E g , decreased by 
increasing the nitrogen concentration, and at a critical 
concentration of cn = 0.35 it disappeared. Fig|3] il- 
lustrates the effects of the nitrogen doping on the E g . 
To clarify our results, we compare the average density 
of states for low and critical nitrogen concentrations. 
FiglU compares the average density of states for the low, 
cn = 0.00005, and critical, cn = 0.35, nitrogen concen- 
tration. At cn = 0.00005, the gap is located at the top 
of the Fermi energy, but inside the pure SWCNT con- 
duction band. By increasing the nitrogen concentration, 
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FIG. 3: The figure shows the effects of nitrogen doping on 
semiconducting gap, E g , for a (10,0) zigzag SWCNT alloy. 
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FIG. 4: The Effects of nitrogen doping on a (10,0) zigzag 
SWCNT,s average density of states. At a fixed boron concen- 
tration cb — 0.15, the average density of states for two nitro- 
gen concentrations, low cn = 0.0005 and critical cn = 0.35, 
are compared. At the critical concentration E g is zero and 
the van How singularities disappear. 

the low edge of the conduction band is moved until the 
gap is closed, hence a semiconductor to semi-metal phase 
transition takes place. Also ID van How singularities dis- 
appeared. 

In the second case, we fixed the nitrogen concentra- 
tion at cn = 0.1, while varying the boron concentra- 
tion. We found that E g , decreases with an increase in 
the boron concentration, and at a critical concentration; 
cb = 0.35005 it tended to zero. Fig \5\ compares the av- 
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FIG. 5: Figure shows the effects of boron doping on the semi- 
conducting gap, Eg, in a (10,0) zigzag SWCNT alloy. 
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concentration, by increasing the nitrogen (boron) concen- 
trations, the Eg decreases and at a critical concentration 
it becomes closed. Therefore, a semiconductor to a semi- 
metal phase transition takes place. Our results show that 
we can control the E g by changing the nitrogen (boron) 
concentration. 

I would like to thanks professor Rafii-Tabar for helpful 
discussion. 
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FIG. 6: Effects of boron doping on a (10, 0) zigzag SWCNT's 
average density of states. At a fixed nitrogen concentration 
of cn = 0.1, the average density of states for two boron con- 
centrations cb = 0.0005 and cb — 0.35005 are compared. At 
a critical concentration, E g is closed and the van How sin- 
gularities disappears and a semiconductor semi-metal phase 
transition takes place. 

erage density of states for the low, cb = 0.00005, and 
critical, cb — 0.35005, boron concentrations. For this 
case, the E g is closed, similar to the first case, and the 
van How singularities also disappeared. Furthermore, the 
semiconductor to semi-metal phase transition was also 
observed. 

In conclusion, we have applied the EMSCA method to 
a disordered SWCNT in order to investigate and role of 
disorder in such materials. For a (10, 10) armchair tube 
and also a zigzag (10, 0) tube we found that the vacancies 
induce some bound states on their host neighbor sites, 
creating a band around the Fermi energy in the average 
density of states. The consequences of this band forma- 
tion around the Fermi energy and also disappearance of 
the ID van How singularities at high vacancy concentra- 
tions is; that the density of states of an armchair and 
also a zigzag SWCNT become similar to a disordered 
(vacancy disorder) 2D graphene sheet density of states. 
A (10,0) zigzag B cb N cn C^ cb ^ cn SWCNT alloy was in- 
vestigated. We found that for a fixed boron (nitrogen) 
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